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Resonant Andreev reflection through superconductor- 
carbon-nanotube devices was investigated theoretically with 
a focus on the superconducting proximity effect. Consistent 
with a recent experiment, we find that for high transparency 
devices on-resonance, the Andreev current is characterized 
by a large value and a resistance dip; low-transparency off- 
resonance devices give the opposite result. We also give evi- 
dence that the observed low-temperature transport anomaly 
may be a natural result of Andreev reflection process. 
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The field of carbon nanotube research has recently en- 
tered a new phase with the fabrication of hybrid device 
structures, in which carbon iiaaotubes are contacted elec- 
trically with other materials^ This is a crucial step, as 
a carbon nanotube-based microelectronics is only possi- 
ble when finite nanotubes can be efficiently fabricated 
and coupled to external leads. In this Letter we re- 
port a theoretical analysis of a hybrid superconductor- 
nanotube junction which has been the subject of a re- 
cent experimental studyu. The experimental device con- 
sists of a single-wall metallic carbon nanotube (SWNT) 
bridging two superconductiM; electrodes. By tuning the 
transpareiicy of the devices, clear signals of Andreev 
reflectionsO were detected via changes in the sub-gap 
resistance at a temperature of T = 4.2K, while other 
transport anomalies were observed at lower T. To date, 
there have been many theoretical and experimental stud- 
ies of normal metal,(M) /superconductor (S) interfaces on 
a mesoscopic scaleETt^. However, no such analysis exists 
for molecular devices. By combining standard-tionequi- 
librium Green's function techniques (NEGF)Bild with a 
tight-binding model (TB) for the SWNT, we have an- 
alyzed quantum transport properties of SWNT-S junc- 
tions. Our results are consistent with the experimental 
datao. 

Although the experimental device consisted of two 
SWNT-S junctions, the data indicates that each of these 
junctions acts independentlyQ. Hence, we focus here 
on the somewhat simpler problem of a N-SWNT-S sys- 
tem, leaving an analysis of the multiple Andreev reflec- 
tions of a S-SWNT-S system for the future. Our theory 
proceeds by combiningJhe NEGF with a standard TB 
model for the SWNTL] such that the coupling of the 
nanotube to the N (left) and S (right) leads are included 
via their appropriate self-energies. By iterating the equa- 
tion of motionEa, standard but tedious algebraO shows 
that there are two contributions to the electric current 



flowing through the device: i.e., I — Ia + Ii- The An- 
dreev current IA^ or sub-gap contribution, is given by 
(^ = e= 1): 

Ia = ^ I dE[fL{E~V)~ fL{E + V)]TA{E), (1) 

where /l,_r denote the Fermi functions of the left and 
right leads respectively, E the electron energy, and V the 
bias potential. The Andreev reflection probability Ta(E) 
is given by: 



Ta{E)= Tr\rLGi2{E)rLG\^{E) 



(2) 



where Tl,r are the appropriate line- width functions de- 
scribing the coupling of the SWNT to the respective 
leads. Ifeua. Gn and G12 are the retarded Green's 
functionsEHrlij of the SWNT which include the proper 
self energies of the leads. These are evaluated by di- 
rect matrix inversioncJ. The remaining contribution to 
the current is given by: 

I,^^ JdE Tr [Ti{E) + T2{E) + T^[E)] 

x[fL{E-V)-fR{E)]pR{E) (3) 

where the density of states of the S lead is pr{E) = 
\E\I^E'^ — for \E\ > A, and zero otherwise; and A is 
the gap energy of the superconductor. Here Ti, T2, T3 
are transmission probabilities for different physical pro- 
cesses that are only nonzero when \E\ > A. Hence, 
these processes describe excitations of the system. In 
particular, Ti = TlGhTrGI-^ gives the familiar tun- 
nefing current; T2 — TLGi2rfjG\2 denotes the branch 
crossing process of Blonder-Tinkham-Klapwijk theoryo; 
and T3 = -{A/\E\)[TlGiiTrG\^ + H.C.] describes the 
Cooper pair formation (or annihilation) inside the su- 
perconducting lead by an incoming electron (or hole). 
Clearly, at zero temperature when bias voltage \ V\ > A, 
all the processes will contribute to current; while for 
|y| < A only the Andreev current I a is nonzero. 

In our numerical calculations, we flxed the gap en- 
ergy at A = 1.45 meV corresponding to that of Nb 
leads; the SWNT TB Hamiltonian is taken to be a 
nearest-neighbor tt— orbital model with bond potential 
Vppjr = —2.75 eV, which is known to give a reasonable 
description of thep£lec.tronic and transport properties of 
carbon nanotubesliZl'Eil. We focus metallic nanotubes of 
flnitc length L. The proximity of the nanotube to the 
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superconductor is modeled through the couphng param- 
eters Tl,r, which are treated as input parameter£3. 

The sohd hues of Fig. (1) show Andreev reflection 
coefiicierit T^C-E) for a metaUic (5,5) SWNT device of 
lengthN L = 3 X 6 + 1 = 19, and for a (9,0) zigzag 
SWNT device of L = 49, with couphng parameter 
Tl ^ Tn ^ 5 X 10~3a.u. (1 a.u. energy is 13.6eV). 
It is clear that resonant transmission with Ta{E) = 2.0, 
i.e., high device transparency, dominates the transport at 
E = 0. This may be understood as follows. It is well- 
known, that infinitely-long armchair nanotubes have two 
states crossing the Fermi level at Ep = 0. On the other 
hand a finite-length isolated armchair nanotube has a gap 
between the two eigenstates near Ep, and this gap isjaaui- 
imized for tube length L = 3n+l where n an integeroEj. 
Therefore, when coupled to the device electrodes, which 
also adds a finite width to the levels, these nanotubes 
have two scattering states at Ep giving rise to large 
transmission with Ta = 2 as shown in Fig.l. The large 
Ta for the zigzag tube is also due to a resonance trans- 
mission through the nanotube states at the Fermi level. 
The two other curves in Fig.l show results for asymmet- 
ric couplings and will be discussed later. Thus, for a 
N-SWNT-S system, the device transparency is critically 
determined by quantum resonance phenomenon which is 
qualitatively different from the case of an infinite SWNT. 
As expected, high transparency devices have a larger Ta, 
since more electrons arriving at the SWNT-S interface in- 
crease the magnitude of the Andreev reflection process. 

Figs. 2 and 3 show the current-voltage characteristics, 
and the differential resistance dV/dl for high- and low- 
transparency (5,5) SWNT devices, respectively. The 
I-V curves of Fig. 2 are qualitatively consistent jjvith 
and quantitatively close to the experimental datao. It 
is clear that a higher slope is observed for the I-V 
curves within the sub-gap range for on-resonance devices 
(dashed curve). This gives rise to an asymmetric resis- 
tance dip near bias voltage V = Q, (see Fig. 3a) with the 
asymmetry due to the finite gate voltage Vg. This resis- 
tance dip is simply a refiection of the high value of Ta for 
on-resonance devices (solid line - Fig.l). The resistance 
dip has a value close to /i/(2 x 4e^) = 2>.2KVi which is pre- 
cisely the expected value of Andreev refiection processes 
in a SWNT-S junciion with two transmitting modes. The 
experimental datalj for two SWNT-S junctions connected 
in series actually gives a value close to b.TKVl, which is 
not far from the expected value of Q.AKQ.. This difference 
is perhaps due to parallel connection of three SWNTs 
bridging i.he superconductor electrodes in the experimen- 
tal setupQ. When the device transparency is low, the dif- 
ferential resistance displays a large peak at y = 0, as 
shown in Fig. 3b. This is consistent both with the I-V 
curve of the low-transparency device shown in Fig. 2 
and the experimental datao. 

So far, the data presented have been for temperatures 
of T = 4.2K, so that features refiecting smaller energy 
scales are completely washed out. However, at^a lower 
temperature of T = 2K, the experimental datao shows 



that a narrow peak emerges in the dV/dl curves for 
zero bias, which is superimposed on the Andreev dip. 
Such anomalous behavior has previously been ascribed 
to the strong electron-electron interactions characteris- 
tic of Luttinger liquids. Surprisingly, our analysis shows 
that these features emerge naturally at lower tempera- 
tures, even within the context of a single-electron theory 
as presented here. This is shown in Fig. 4, which illus- 
trates the emergence of a narrow peak out of the overall 
Andreev dip as the temperature is lowered. 

We can trace this low temperature anomaly to the 
basic physical process which gives rise to the Andreev 
current. To demonstrate this, we neglect complications 
due to the molecular structure of the nanotubes and 
assume that the resonant Andreev process is mediated 
by a single state at energy Eq < A, i.e. we "shrink" 
the nanotube to a simple quantum well with a sin- 
gle level. This is qualitatively reasonable because the 
sub-gap energy scale is set by gap energy A which is 
much-smaller than the level spacing of the nanotube we 
studycJ. Hence, we expect that only the two degener- 
ate levels at the Fermi energy will contribute apprecia- 
bly to the Andreev current. For this single- level case, 
the Green's functions are drastically simplified^ so that 
equivaleiitly, familiar scattering matrix theory can be 
appliedEl. Near the resonance, the transmission ampli- 
tude in the normal state assumes the Breit-Wigner form, 
t{E) = iy/TLTRj{E -Eo + ir/2); and the reflection am- 
plitude becomes r{E) = 1~ iTii/{E — Eo + «r/2); where 
r = r^+r PI is the total line width. For these amplitudes, 
Eq. (0) takes on the Breit-Wigner form: 



Ta{E) 



"p2 T^2 



4[ij2_^2+r|r] 



^2 [r2 + 6T^] ' 



(4) 



where ST = Tp—Tfj. This result allows us to draw several 
conclusions. For simplicity, we set the level at Eq = 0, to 
obtain 



TAiE) 



A{E^ 



^5^ ^ 

4 J 



(5) 



These expressions reduce to tho-result of scattering ma- 
trix theory for N-Dot-S systemsO when we set E = Eg = 
0. Eq.(|) indicates that if Vp = Tr so that (JF = 0, 
then resonant Andreev reflection occurs at E = with 
Ta{E = 0) ^ 1. On the other hand, if > F^^ so 
that ST > 0, Ta takes on a maximum value at i? = 
but this maximum value is less than 1. For nanotubes 
this situation is shown by the dotted lines of Fig.l. Fur- 
thermore, if Fl < Tp such that JF < 0, Ta is charac- 
terized by two resonant peaks with Ta = 1 at energies 
E± — ±\/—rST/2. For nanotubes this behavior is shown 
by the dashed lines of Fig.l. Hence, due to a split be- 
tween the electron and hole levels when the nanotube 
is in contact with a superconducting lead, the Andreev 
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coefficient in the sub-gap region can display different be- 
haviors. When Ta{E) displays a double peak, the An- 
dreev current Ia shows a resistance anomaly such that 
a small peak develops inside the overall dip at low tem- 
peratures. At higher temperatures such as A.2K, the 
anomaly is smeared out, and hence, not observable. One 
can also confirm that, qualitatively, the above conclusion 
holds for cases of nonzero gate voltage. 

The analysis presented so far has been for (5,5) nan- 
otubes witii_a length oiL = 19 unit cells, and (9,0) zigzag 
nanotubesE3 with L = 49. However, we expect our results 
to be general for other mctallic-|nanotubes and lengths. 
Previous investigations showecO that transport through 
finite armchair nanotubes differ qualitatively from the 
infinite length limit. In particular, if L = 3n + 1 where 
n is an integer, the armchair tube has large conductance 
due to the crossing of scattering state energy levels at the 
Fermi energy. Other tube lengths produce much smaller 
conductance due to a gap between the scattering states. 
The inset of Fig.2 shows the I-V curve of N-SWNT-S 
devices with L = 17 and L = 18: these devices have 
very small currents because their device transparency 
are drastically diminished by the energy gap between the 
scattering states of the SWNT. Experimentally, such an 
energy gap along with conductance oscillations on finite 
tubes have already been detected with scanning probesEj. 

In summary, we have investigated the sub-gap trans- 
port properties of N-SWNT-S systems and our results are 
consistent with a number of experimental observations. 
Specifically, the dependence of the Andreev current on 
the device transparency, the behavior of the differential 
resistance in the sub-gap region, as well as the observed 
low-temperature resistance anomaly may all be explained 
in terms of the proximity of the nanotubes with the su- 
perconducting lead. There are, however, still several is- 
sues which cannot be studied within the context of our 
model. First, the experimental dataci showed a sensi- 
tive gate voltage dependence at very low temperatures 
on the order of 40 mK, where the differential resistance 
anomaly could be a peak or a dip depending on the value 
of the gate voltage. This behavior has been jconsidered 
as likely due to electron-electron interactionsa. Second, 
there are experimental and theoretical evide nqe t hat nan- 
otubes can have non-Fermi liquid behaviornZlna. It will 
be of great interest to investigate the situation of a non- 
Fermi liquid model of nanotubes in contact with a super- 
conductor to see if other, finer features emerge. Another 
important problem is the detailed atomic structural anal- 
ysis of the nanotube-superconductor interface which, to 
a large extent, controls the interface transparency. Fi- 
nally, although we do not expect charge transfer to play 
a critical role in understanding the Andreev current for 
SWNT-S interfaces because of the superconducting gap, 
a more complete investigation of this delicate effect will 
certainly enhance our understanding of quantum trans- 
port for nanoscale devices. 
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FIG. 1. Andreov reflection coefBcient Ta as a function of 
electron energy E. Main graph: for devices consisting of (5,5) 
SWNT with length L=19 unit cells. Inset: For (9,0) zigzag 
nanotube systems. Other parameters are fixed as zero bias 
and zero gate voltages. Solid curve: Tl ~ Pij = 5 x 10~^; 
dotted curve: Fl = 6 x 10"^, Tr = Ax 10"^; dashed curve: 
Fi = 4 X 10"^, Fh = 6 X 10"^. Here the Fs are measured in 
a.u.. 

FIG. 2. LV curves for the N-SWNT-S device at temper- 
ature 4.2K. The SWNT is a (5-5) metallic tube with length 
1/ = 19 unit cells. Solid curve is for a low transparency device 
(off resonance transmission) with parameters Fz, = Fjj = 0.8 
and Vg = 0. Dashed curve is for a high transparency de- 
vice (on resonance transmission) with Fz, = F_r = 5 x 10"'^ 
and Vg = 0.6mV. Inset: I-V curves (same units as the 
main graph) for tubes with L = 17 and 18 unit cells with 
Fz, = Fh = 5 x 10"^ and Vg = 0.6mV. 

FIG. 3. Differential resistance dV/dl as a function of bias 

voltage V. (a) For high transparency (on resonance) device 
corresponding to the dashed I-V curve of Fig. 2. (b) For low 
transparency (off resonance) device corresponding to the solid 
LV curve of Fig.2. (5,5) SWNT with length L = 19 unit cells 
are used. 

FIG. 4. Differential resistance dV/dl for high trans- 
parency devices at different temperatures. A peak emerges 
from the overall dip as temperature is lowered, (a). For (5,5) 
nanotube devices with L = 19 unit cells, (b). For (9,0) 
nanotubes with L = 49. Other parameters: Vg = 0.6mV, 
Fz, = 3 X 10"^, Fz{ = 8 X 10"^. 
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